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A supercurrent switch in graphene 7r -junctions 
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We study the supercurrent in a superconductor/ferromagnet/superconductor graphene junction. In contrast to 
its metallic counterpart, the oscillating critical current in our setup decays only weakly upon increasing exchange 
field and junction width. We find an unusually large residual value of the supercurrent at the oscillatory cusps 
due to a strong deviation from a sinusoidal current-phase relationship. Our findings suggest a very efficient 
device for dissipationless supercurrent switching. 



PACS numbers: 73.23.Ad, 74.50.+r,74.45.+c, 74.78.Na 

Graphene is a condensed matter system displaying an emer- 
gent low-energy 'relativistic' electronic structure. For un- 
doped graphene, the Fermi level reduces to six points, giving 
rise to nodal fermions at the edges of the Brillouin zone. Cur- 
rently, it is of considerable interest, and potentially of great 
technological importance, to investigate how such unusual 
low-energy electronic structures manifest themselves in het- 
erostructures where proximity effects are prominent. In par- 
ticular, potential for future applications in devices seems plau- 
sible if such proximity-structures would combine two major 
functionalities in materials science, namely magnetism and 
superconductivity. 

Recently, there have been several experimental reports on 
proximity-induced superconductivity in graphene yj, |2||. A 
measurable supercurrent was observed between regions of 
graphene under the influence of proximity-induced supercon- 
ductivity in these works. Due to the massless nature and 
energy-independent velocity of the charge-carriers, graphene 
offers a unique environment for the manifestation of a Joseph- 
son effect. Unusual behavior for the supercurrent in a super- 
conductor/normal/superconductor (S/N/S) graphene setup has 
been predicted, including an anomalous scaling behavior with 
the length of the normal region in the undoped case [3] and an 
oscillatory behavior as a function of gate voltage in the normal 
region fl]. 

Interestingly, it has also been shown ||5l |6[] that ferromag- 
netic correlations may be induced in graphene nanoribbons 
by means of external electrical fields. A suggestion for a more 
conventional magnetic proximity effect has also been put forth 
S^ i 03,181], by means of exploiting a magnetic gate in contact with 
a graphene layer. The accompanying exchange splitting be- 
tween the spin-f and spin-], electrons in graphene has been 
estimated B7[] to lie around 5 meV for the magnetic insula- 
tor EuO. Precise estimates of the proximity induced exchange 
splitting are difficult in this case, due to the strong effect of 
the proximity layer on the magnetization in EuO |9|. Nev- 
ertheless, it is known that the magnetization in the proximity 
EuO layer is tunable [|lfjf] . In applications, this is a great ad- 
vantage, since it in principle offers the possibility of a tunable 
proximity induced magnetization in graphene. Moreover, re- 
cent experiments on spin injection in a graphene layer show 
a rather long spin relaxation length ~ 1 ^m at room tempera- 
ture. This indicates that graphene is a promising material for 
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a) The Andreev bound state e+ [Eq. l[3}] for A0 = 7r/2. 
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b) The Andreev bound state e_ [Eq. <[3j] for A(f> — tt/2. 

FIG. 1 : (Color online) Contour-plot of the Andreev bound states in 
the ferromagnetic region carrying the current between the supercon- 
ductors. 

The following question arises naturally: Do novel physi- 
cal effects arise due to the peculiar electronic properties of 
graphene and simultaneously the interplay between ferromag- 
netic and superconducting correlations? The wide range of 
exotic phenomena that originate from the mutual interplay be- 
tween magnetic and superconducting order include Q-tt tran- 
sitions I12II13I1 . odd-frequency pairing B14I1 . and even intrin- 
sic coexistence of ferromagnetism and superconductivity in 
the same material J15[|l6ll . In particular, from the viewpoint 
of applications, the possibility of altering the fundamental 
Josephson current-phase relationship in a controlled fashion 
may bring about potential implications for their use in super- 
conducting electronics as well as in (quantum) logic circuits 
based on superconductors B17I1 . 

In this Letter, we investigate the interplay between 
proximity-induced superconductivity and ferromagnetism in 
a graphene layer, resulting in an unusual behavior of the su- 
percurrent through the system. Our main results are: i) The 
current-phase relationship deviates strongly from sinusoidal 



behavior, indicating a significant contribution from higher har- 
monics and ii) the critical current at the 0-ir transition is finite 
and has a much larger value than the one observed in metal- 
lic systems. The latter result suggests a very efficient perfor- 
mance of the device as a supercurrent switch. 

We envisage an experimental setup where superconductiv- 
ity is induced in two parts of the graphene region by means 
of conventional superconductors, such as Nb or Al, in close 
proximity. Between the superconducting regions, an exchange 
splitting is induced in the graphene layer by means of e.g. a 
magnetic insulating material. Instead of using a magnetic in- 
sulator such as EuO, where one in principle could tune the 
magnetization in the proximity layer with an external mag- 
netic field, one also could envision using a multiferroic (e.g. 
BiFeOa) or piezomagnetic material (e.g. Fe a; Nij / B 2 ) in close 
proximity to the graphene layer. Both of these classes of mate- 
rials would offer the opportunity of tuning the exchange field 
in the material by some external control parameter - electric 
field due to the magnetoelectric coupling in the former case, 
and pressure in the latter. Upon modifying the exchange field 
in the proximity layer of the material, it is reasonable to as- 
sume that the proximity-induced exchange field in graphene 
would also be altered. Materials in which the magnetoelectric 
coupling is substantial are currently attracting much interest 
due to their potential for novel technological applications [18]. 
In order to control the local Fermi level in the ferromagnetic 
(F) region, one could possibly use a normal gate on top of the 
magnetic insulator to create a tunable barrier [7]. The super- 
conducting (S) regions are assumed to be heavily doped, such 
that the Fermi energy satisfies ep ^> A, while the F region is 
taken be undoped, i.e. e' F ~ 0. Moreover, we assume sharp 
edges for the region separating the F and S graphene regions, 
and focus on the short-junction regime which is experimen- 
tally feasible. 

We will proceed to show that the Josephson current in an 
S/F/S graphene junction displays a strong oscillatory, non- 
monotonic dependence on both the exchange field h and width 
d of the junction. Most interestingly, we find a large residual 
value of the supercurrent at the cusps of these oscillations. 
This indicates a sign reversal of the current, and the consider- 
able residual value of the supercurrent at these cusps suggests 
a very efficient device for dissipationless supercurrent switch- 
ing. We now present our results in detail. 

The F region separating the superconductors is taken to be 
undoped, such that the effective Fermi level is ah for spin-er 
electrons. The regions S must be strongly doped to justify the 
mean-field treatment of superconductivity. We assume that 
this is comparable to the estimated exchange-splitting in the 
F region y[ |8|]. Thus, we take ep — h to obtain analyti- 
cally tractable results. To construct the scattering states that 
carry the supercurrent across the F region, we write down the 
Bogoliubov-de Gennes equations |19] in the presence of an 
exchange field h. The Bogoliubov-de Gennes equation essen- 
tially describes the eigenstates of quasiparticles in each of the 
graphene regions and their belonging eigenvalues e. It may 
be obtained by diagonalizing the full Hamiltonian, and con- 



stitutes the foundation for constructing the scattering states 
which are involved in the transport formalism we here use. 
For the spin-species a, one finds that 



Hq — ah{x) <jA{x) 
crA*(x) -H -ah(x] 
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Here, we have made use of the valley degeneracy and de- 
fined Hq = vpp • c, where p is the momentum vector in 
the graphene plane and er is a vector of Pauli matrices. The 
superconducting order parameter A(x) couples electron- and 
hole-excitations in the two valleys located at the two inequiv- 
alent corners of the hexagonal Brillouin zone. The u a spinor 
describes the electron-like part of the total wavefunction ip a = 
(w <T , w _<T ) T , and in this case reads u a — {^f A + ,ip% + ) T while 
v~ a = Tu a . Here, T denotes the transpose while T is the 
time-reversal operator. To capture the essential physics, we 
write A (a;) = Aoe 1 * LR in the left and right S region and 
A(x) = otherwise. Similarly, we set h(x) — h in the F re- 
gion and h = otherwise. The Josephson current is computed 
via the usual energy-current relation summed over projections 
of all paths perpendicular to the tunneling barrier [20] 
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where £i(A(j>) are the Andreev bound states carrying the cur- 
rent in the F region, and A0 = 0r — </>l is the macro- 
scopic phase difference between the superconductors. The 
integration over angles 7 takes into account all possible tra- 
jectories and f(x) is the Fermi-Dirac distribution function. 
We define the critical supercurrent as I c = \max{Ij(A<fi)}\ 
and introduce Iq = 2eAo. The procedure for obtaining 
Ei(A<p) is the same as in Ref. [3] and the details will be 
given elsewhere; here we give the main results. By in- 
troducing T(-/,P) = 4 sin 2 7 sin 2 (P cos 7) + cos 4 7 and 
$(7,P) = [2 sin 2 (P cos 7) — cos 2 7] cos 2 7COS A0— cos 4 7— 
4 sin 2 7sin 2 (Pcos7), we find that the allowed bound states 
have energies ±e cr (A0) (er = ±) with 
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The parameter P = hd/vp captures the effect of both the ex- 
change field h and the length d of the junction. To understand 
the nature of these bound states, consider Fig. Q]for a represen- 
tative plot of £±(A</>), using A<fi — tt/2. As is seen from both 
plots, the bound state energies exhibit a strong oscillatory de- 
pendence on the parameter P. This indicates that similar os- 
cillatory behavior may be expected in the supercurrent itself. 
Interestingly, the oscillations seen in Fig. Q]are not damped 
with increasing P. This directly reflects the Dirac-cone linear 
dispersion of the graphene electrons and is reminiscent of the 
weak damping of conductance oscillations at subgap energies 
in graphene-superconductor junctions 112 ll 12211 . 



Inserting the derivative of Eq. (01 into Eq. (fJJ provides 
the supercurrent. The current-phase relationship for the S/F/S 
graphene junction is shown in Fig. [2] With increasing P, the 
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FIG. 2: (Color online) a) Current-phase relationship in the S/F/S 
graphene junction with undoped F region. We have fixed ep = h 
and set e F — 0. We have used values of P in the interval [0.0, 1.2] 
in steps of 0.2. b) Current-phase relationship in the S/F/S graphene 
junction with doped F region. We have set £f/Ao = 10, e F /Ao = 
15, d/£ = 0.05 and vary h/Ao in the range [0, 5] in steps of 1. 

critical current gets suppressed and finally the sign of the cur- 
rent is changed. Remarkably, the critical current never goes 
to zero. An interesting feature of the plot in Fig. |2^) is 
that the discontinuity at A<j) = ir for P = is split for in- 
creasing P. The discontinuity of the current-phase relation 
originates with a crossing of Andreev levels in the normal 
graphene (F graphene with h = 0) region at A<f> = tv. For 
Acf> £ [0, 7r), only the 0-mode Andreev bound state carries the 
current. For A<j> £ (it, 2ir], the 7r-mode Andreev bound state 
carries the current, such that there is an abrupt crossover ex- 
actly at Acf> = ir. The situation changes when h ^ 0, since the 
spin-splitting doubles the number of Andreev bound states. 
Consequently, the crossover between different modes may oc- 
cur at Acj) ^ 7T, as a result of the superharmonic current-phase 
relationship. We have checked explicitly that the strong de- 
viation from a sinusoidal current-phase relationship persists 
for larger d that do not satisfy d/£ <C 1. However, in this 
case one should strictly speaking also include the contribu- 
tion to the current from the continuum of supergap states J20I1 . 
This requires a separate study, and we here focus on the short- 
junction regime. 

To show that the splitting of this discontinuity originates 
with the presence of an exchange field which separates the 
spin-f and spin-j bands, we have also numerically solved the 
current-phase relationship for a nonzero Fermi level in the 
ferromagnetic region. Although we have obtained analytical 
results in this regime, these are somewhat cumbersome and 
therefore omitted here. The result is shown in Fig. |2j)) where 
we have chosen Ao = 1 meV, ep = 10 meV, and e' ¥ = 15 
meV, and varying h in the range [0, 5] meV. This ensures that 
there are no evanescent modes, such that only the Andreev 
bound states carry the current. We choose a junction with 



d/£ = 0.05, where £ is the superconducting coherence length, 
since the short junction regime d <C £ is the experimentally 
most relevant one. The figures in a) and b) correspond to two 
quite different regimes: in a) the exchange field is much larger 
than the Fermi level while in b) the exchange field is much 
smaller than the Fermi level. The trend is nevertheless seen 
to be the same in both cases, namely a progressive splitting of 
the discontinuity located at Ac/) = it in the paramagnetic case. 
Assuming a heavily doped superconducting region with 
£f = 10 meV and an effective gap Ao of 1 meV, a mean- 
field treatment is justified by e-p ^S> Ao. Moreover, the short- 
junction regime requires that d <C £■ Using vp ~ 10 6 m/s 
in graphene, we obtain from £ = vp/A that d <C 650 nm 
is required. This condition has been met in at least two ex- 
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FIG. 3: (Color online) The critical supercurrent in a proximity- 
induced S/F/S graphene junction for sf = h and e F = 0. 

perimental studies of proximity-induced superconductivity in 
graphene fljj, |2J]. The critical supercurrent I c for an S/F/S 
graphene junction for ep = h and e F = is shown in Fig. 
[3] The critical current shows oscillations with respect to P, 
but decays weakly compared to the metallic case and never 
reaches I c = in the relevant regime. For instance, there is a 
factor ~ 100 in reduction of the amplitude of the current right 
after the second cusp in the metallic case for h ~ 10Ao (see 
Fig. 2 in Ref. 23) while there is only a factor ~ 2 in reduction 
of the amplitude in the present case. Right at the cusps lo- 
cated at P ~ {0.8, 2.8, 4.4}, there is a large residual value of 
the supercurrent which should be experimentally detectable. 
This is very distinct from the usual sinusoidal current-phase 
relationship for the Josephson current, in which the supercur- 
rent vanishes completely at the 0-7T transition. The first switch 
occurs at a value P = hd/vp ~ 0.8. For an exchange split- 
ting of h ~ 10 meV, this requires a junction width d = 50 
nm. Alternatively, employing a junction width of d = 100 nm 
1 1 , \M one would need an exchange splitting of h ~ 5 me V 
0. 

In order to explain the appearance of cusps in the critical 
current dependent on exchange field and junction width, it is 
instructive to draw parallels to the metallic S/F/S junction and 
the behavior of the supercurrent. In most experimental sit- 
uations, the effective barriers separating the F and S regions 



are strong, leading to a current-phase relationship which is 
very nearly sinusoidal, i.e. I c = IosmAtfi [24]. By tuning 
the temperature T and width of the junction d, one is able to 
switch the sign of the amplitude I , which necessarily means 
that one must have 1$ = at some point. Precisely such 
behavior has been observed in several experiments 1U2U25I1 . 
In the present system, the current-phase relationship deviates 
strongly from sinusoidal behavior, and contains a significant 
contribution from higher harmonics. Tracking the absolute 
value of the current with increasing P from Fig. it is seen 
that I c never becomes zero. Instead, it has a large residual 
value at the points where the current changes sign. While a 
small, but finite value of the supercurrent at the 0-ir transi- 
tion also has been observed in metallic S/F/S junctions M26I1 . 
the magnitude of the residual value of the supercurrent in the 
graphene case is huge compared to the metallic case. 

Since we have assumed a homogeneous chemical poten- 
tial in each of the S and F graphene regions, the experimen- 
tal observation of the predicted effects require charge homo- 
geneity of the graphene samples. This is a challenge, since 
electron-hole puddles in graphene imaged by a scanning sin- 
gle electron transistor device [27] suggest that such charge in- 
homogeneities play an important role in limiting the transport 
characteristics of graphene close to the Dirac point B28I1 . In 
doped graphene, as considered here, we expect that the in- 
homogeneities should play a smaller role than in undoped 
graphene. Although we have neglected the spatial variation 
of the superconducting gap near the S/F interfaces, we do not 
expect our qualitative results to be affected by taking into ac- 
count the reduction of the gap. Also, we have assumed that 
the junction d is short enough to neglect the orbital effect the 
magnetic field constitutes on the electrons. 

In summary, we have investigated the interplay between 
proximity-induced superconductivity and ferromagnetism in 
a graphene layer. In contrast to its metallic counterpart, the 
oscillating supercurrent in our setup decays only weakly upon 
increasing exchange field and junction width. We find huge 
residual values of the supercurrent at the 0— 7r transition points 
where the supercurrent changes sign. If the exchange splitting 
could be adjusted by means of some external source, such as 
gate voltage or external electrical fields yl |6j, ll8[], these re- 
sults imply that the supercurrent across the junction may be 
tuned in a controllable fashion. Specifically, a very efficient 



supercurrent switch may be realized by tuning the exchange 
field infinitesimally near the sign reversal points. 
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